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Abstract 

We analyze dynamics of a quantum particle in a square lattice in the Hall configuration beyond 
the single-band approximation. For vanishing gauge (magnetic) field this dynamics is defined by 
the inter-band Landau-Zener tunneling, which is responsible for the phenomenon known as the 
electric breakdown. We show that in the presence of a gauge field this phenomenon is absent, at 
least, in its common sense. Instead, the Landau-Zener tunneling leads to appearance of a finite 
current which fiows in the direction orthogonal to the vector of a potential (electric) field. 
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I. INTRODUCTION 



This work continues our studies of the Landau-Zener tunnehng (LZ-tunnehng) for a 

n 

quantum particle in the Hall configuration To be certain, here we assume a charged 
particle in electric and magnetic fields, although the results are equally applied to a neutral 
particle (for example, an atom in an optical lattice) subject to artificial gauge and potential 
fields. Specifically, we consider the following system 
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H 



Px + iPy - -BxY 



+ V{x,y) + eiF^x + Fyy) , (1) 
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where B is the magnitude of a magnetic field, F = {F^^Fy) an electric field, and V{x,y) = 
Vx cos(27rx/a) +Vy cos(27r?//a) the periodic potential. For 5 = the spectrum of the system 
(II]) consists of Bloch bands separated by energy gaps. It is well know that a strong electric 
field induces Landau-Zener transitions (tunneling) across the energy gaps. In our previous 
work pj] we have shown that these transitions also take place in the presence of a magnetic 
field and calculated the tunneling rate. Unfortunately, the analytical approach of Ref. |l| 
provides only the rate of tunneling but is incapable to describe dynamics of the system ([T]) 
if LZ-tunneling is not negligible. In the present work we study dynamical manifestations 
of LZ-tunneling by numerically solving the Schrodinger equation with the Hamiltonian ([1]) 
and comparing the results with those obtained on the basis of the tight-binding model (tb- 
model), 

(^tfe^)^ ^ = -y (^z+i,™ + ^z_i,™) - ^ (7/>i,^+ie^2™z ^ ,pi^^_^e-'^--') +ea{FJ + Fym)iji,^ , 

(2) 

which is believed to approximate the Hamiltonian ([1]) in the case of negligible LZ-tunneling. 
In Eq. ([2]) J^^y are the hopping matrix elements in two orthogonal directions, a is the lattice 
period, I = x/a and m = y/a label lattice sites, and a is the Peierls phase, 

" = 1^- 

It should be mentioned that the tight-binding approximation implicitly assumes a deep 
lattice. For the considered periodic potential this requires Vx,y > Er/2, where Er = h^/Ma^. 
If the later condition is satisfied the ground Bloch band of the one-dimentional Hamiltonians 

ffr^^-«.-(2.£). (4) 




FIG. 1: Squared absolute value of the continuous wave-function ^{x, y, t), left, and squared absolute 
values of the amplitudes ipi^mit), right, for t = 300. The system parameters are Vx = Vy = 0.5, 
F = 0, and B/Itt = a = 1/8. 

{Hq^^ is given by the similar equation) are well approximated by the cosine function, 
E''^\kx) = Eq^^ — JxCos{aK,x)- Thus the next to neighboring hopping terms can be in- 
deed neglected in the Hamiltonian (|2]). In what follows we measure the energy in units of 
Er, the length in units of a (which we set to 2tt), and drop all physical constants. 



II. VANISHING ELECTRIC FIELD 

We begin with comparing dynamics of the systems ([1]) and ([2]) for F = 0. We solve 
the time-dependent Schrodinger equation with the Hamiltonian ([1]) on the grid and check 
the convergence of the results with respect to the space and time discretization. Having 
the wave-function obtained, we calculate its overlap with the Wannier functions ^i^m{x,y) 
associated with the ground Bloch band of the Hamiltonian Hq = Hq^^ + Hq^^: 

^i,m{t) = J dxdy^i^rn{x, y)'^{x, v, t) , ^i,m{x, v) = ^oflix - 2itI, y - 2nm) . (5) 

Finally the probability amplitudes are compared with those calculated on the basis of 
the tight-binding Hamiltonian ([2]). 

A remark concerning initial conditions is in order. To facilitate comparison with tb- 
model we choose \l/(a;, = 0) = $o,o(a;,2/), that corresponds to population of the single 
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FIG. 2: Upper panel: participation ratio ([6]) as functions of time for the parameters of Fig. [TJ 
Lower panel: the same yet F = 0.015, F^/Fy = - l)/A ^ 0.309. Dashed lines show the 
participation ratio calculated by using tb-model. 

site, ipi^mit = 0) = 6ifi6m,o- Notice that this initial state assumes population of the whole 
ground Bloch band if S = or all ground magnetic bands ii B ^ l|]. As an example 
the left panel in Fig. [1] shows the squared absolute value of the continuous wave-function 
\E'(x, y, t) for the specified initial condition after evolution time of the order of one cyclotron 
period T^. = 27r /ujc [2^. The right panel in Fig. [1] shows the squared absolute values of the 
amplitudes ([5]). For the considered evolution time the depicted in the right panel populations 
of the lattice sites are undistinguishable by eye from those obtained on the basis of tb-model. 
However, for larger times we observe clear deviations between the original system and tb- 
model, which we discuss in the next paragraph. 

As an overall characteristic of the system dynamics we consider the participation ratio 

Pit) = (T.\AUt)A . (6) 

The dashed hne in the upper panel in Fig. [2] depicts P{f) calculated by using tb-model. The 
participation ratio shows oscillatory dynamics, superimposed with a slow increase in its mean 
value. This slow increase is due to finite widths of the magnetic bands, which imphes ballistic 
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FIG. 3: The projection ([7]) on the subspace spanned by the ground WS-states for F = (sohd 
hne) and F = 0.015 (dashed hne). The dash-dotted hne shows the decrease of total probabihty in 
the case F = 0.015 due to absorbing boundary conditions. 

dynamics for asymptotic times ^|. The dashed line should be compared with the solid line, 
which is calculated by using the original Hamiltonian. Although qualitative correspondence 
is noticed, quantitative agreement holds only for finite time. The reason for the observed 
deviations is a relatively large value of a used in our numerical simulations. In fact, if i? 7^ 
the condition v^^y > 1/2 alone is not enough to justify the tight-binding Hamiltonian. The 
other necessary condition, which is often overlooked when deriving ([2]) from the continuous 
Hamiltonian ([1]), is a <^ 1. Indeed, if a 7^ dynamics of the system ^ is not restricted 
to the subspace of the Hilbert space spanned by the ground Wannier states y), even 

if the initial wave-function belongs to this subspace. This statement is illustrated in Fig. [3] 
which shows the projection of the wave function on the specified subspace, 

N{t) = Y.\i^,At)\' . (7) 

Obviously, the smaller a = 27tB, the closer N{t) is to unity. However, smaller a require 
larger lattices and longer evolution time to see characteristic features of the system dynamics. 
Our choice a = 1/8 is a compromise between accuracy of the tight-binding approximation 



and complexity (computational time) of numerical simulations. 



III. VANISHING MAGNETIC FIELD 

Next we consider the case F 7^ and B = 0. For B = dynamics of the system ([1]) are 
Bloch oscillations, where the mean position (normal mode) or the width (breathing mode) 
of a localized wave packet oscillates with the Bloch frequencies 

= eaF,^y/h = 27rF,, . (8) 

Simultaneously, the oscillating wave packet emits sub-packets, which get accelerated by the 
electric field and move away from the main packet along the crystallographic axes of the 
lattice . This dynamics is exemplified in Fig. m which shows |\E'(x, t) p for Fy = 0.015 
and = 0, after evolution time t = 10000. For the chosen electric field configuration the 
sub-packets move in the negative y direction while the packet itself spreads in the x direction. 
Since in numerical simulations we are forced to deal with finite lattices, we impose periodic 
boundary conditions at x = ±16 and absorbing boundary conditions at y = —16. For 
this setup the system quickly reaches quasi stationary state, where time-average occupation 
probabilities of the lattice sites do not depend on /. Notice that due to absorbing boundary 
condition the total probability decay in time. The dash-dotted and dashed lines in Fig. [3] 
show the total probability and the projection ([7]) on the subspace spanned by the ground 
WS-states. These two curves are seen to follow each other. Thus, unlike the case of null 
electric field, deviations from tb-model are now exclusively due to irreversible LZ-tunneling 
to higher bands. 

Along with the quantity (JTj) the other important characteristics of the system is the mean 
kinetic energy Exit) of the particle. Since the total energy is conserved, an increase in the 
kinetic energy is compensated by decrease in the Stark energy 

Es{t) = J dxdy\^{x,y,t)\\F,x + Fyy) . (9) 

It is easy to prove that E^ ~ —Es{t) infinitely grows due to LZ-tunneling. This is in strong 
contrast with the case B ^ 0, which we shall study systematically in the next section. 
For the moment we only mention that a finite B prohibits infinite increase/decrease in the 
kinetic/Stark energy and imposes an upper boundary for possible energies of the quantum 
particle. This boundary can be estimated by using the classical arguments as follows. 
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FIG. 4: Logarithm of y, t)p at t = 10000 for vanishing magnetic field. The system parameters 



are Vx = Vy = 0.5, Fy = 0.015, and = 0. 
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FIG. 5: Classical trajectories of 10 particles with random initial conditions in the central potential. 
Parameters are Vx = Vy = 0.5, B = 0.02, Fy = 0.02, and Fx = 0. 
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Let us consider classical counterpart of the Hamiltonian ([T]) where, for the sake of simplic- 
ity, we assume = 0. Let us also assume that the particle is initially located at {x,y) ~ 
and its energy is slightly higher than the potential barrier in the y direction. Then, for finite 
time defined below, the particle is accelerated by the electric field in the y direction and 
Py(t) ~ Ft. During this initial stage the motion in the x direction is oscillations around the 
central local minimum of the effective potential V{x,t) = —VxCosx + {Ft — BxY /2. It is 
easy to see that this local minimum disappears if t > t*, where 

At this time the particle escapes the central valley of the 2D potential with the kinetic 
energy 

E^a. ~ {Fef/2 ~ {vx/Bf . (11) 

As soon as the particle escapes the valley its further trajectory becomes usual cycloid, see 
Fig. O where the kinetic and Stark energies oscillates around Emax- The above analysis can 
be easily generalized for other directions of the electric field vector, resulting in the same 
estimate. It should be also mentioned that Eq. (fTOl) gives the maximal escape time. If we 
consider an ensemble of particles with slightly different initial conditions, we get different 
escape times and, as the consequence, a broad distribution for the kinetic/Stark energies. 



IV. FINITE MAGNETIC AND ELECTRIC FIELDS 

We proceed with the case of F 7^ and B 0. First we consider the situation where F is 
aligned with the y axis. Figure [H] shows ?/, t)p for F = 0.015 and a = 1/8 (S ~ 0.02) 
and should be compared with Fig.Hlwhere B = 0. Qualitative difference between two figures 
is that in the former case the tunneling particle acquires finite velocity in the x direction, 
which is consistent with the classical dynamics considered in Sec. IIIII We calculated the 
mean current along the x direction, which was found to differ from zero. Notice that tb- 
model predicts zero current for the considered initial condition. Thus a finite current in the 
X direction can be considered a manifestation of LZ-tunneling. 

Unfortunately, we have not been able to find the actual value of the current. To do this 
one should move the absorbing boundary further to negative y till the convergence is reached. 
[According to the estimate (fTOj) this requires the lattice of the order of 1000 sites, which 
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FIG. 6: Logarithm of at t = 10000. Parameters are Vx 



0.5, a = 1/8 {B « 0.02), 



Fy = 0.015, and = 0. 



is beyond our computational capabilities.] In this case we also expect to see a qualitative 
difference in the tails of the reduced probability 



piy) 



|^(x,y,t)pdx 



(12) 



see Fig. [TJ which defines the Stark energy ([9]) trough the equation Es = F J yp{y)dy. Based 
on the classical consideration, if i? 7^ the reduced probability is expected to decay faster 
than l/\y\ which ensures finite kinetic and Stark energies. On the contrary, if B = the 
reduced probability should show a diverging tail because the Stark and kinetic energies 
infinitely increase when t ^ 00. 

Finally we consider the case of an irrational direction of the electric field F^/Fy = {\/5 — 



l)/4 ^ 0.309. As shown in the recent works [7|-^ devoted to dynamical and spectral 
properties of the system ([2]), in this case the electric field localizes the particle in the lattice, 
where the localization length tends to one site for F/J ^ 1. Notice that this result holds 
for both zero and finite magnetic fields. Numerical simulations of the original system for 
F = 0.015 partially confirm this conclusion: When plotting the wave-packet in the linear 
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FIG. 7: Reduced probability p{y) at t = 10000 in the logarithmic scale. The solid and dashed lines 
correspond to a = 1/8 and a = 0, respectively. 

scale we observe that only 3-4 sites participates in the dynamics. Moreover, if we normalize 
the amplitudes ([5]) against N(t) the participation ratio ([6]) practically coincides with that 
calculated on the basis of tb-model, see lower panel in Fig. O However, when plotting the 
wave packet in the logarithmic scale we clearly see the effect of LZ-tunneling, as well as the 
qualitative difference between the cases 5 = and B 0, see Fig. [HI 

V. CONCLUSIONS 

We analyzed LZ-tunneling for a charged particle in a square lattice subjected to normal 
to the lattice plane magnetic field B and in-plane electric field F. For B = the rate of LZ- 
tunneling is defined by the size of the energy gap A separating the ground Bloch band from 
the rest of the spectrum and the magnitude of the electric field. If the rate of this tunneling 
is not negligible (strong electric fields), the ground Bloch band rapidly depletes and the 
particle is accelerated towards infinite kinetic energy - the phenomenon known as electric 
breakdown. We showed that a finite magnetic field essentially modifies this phenomenon, 
where the effect of the magnetic field is twofold. 
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FIG. 8: Logarithm of |^'(x,y,t)p at t = 10000 for a = (left panel) and a = 1/8 (right panel). 
The other parameters are Vx = Vy = 0.5, F = 0.015, and F^/Fy = 

First, the magnetic field splits Bloch bands into magnetic bands, thus introducing the 

n 

new energy gaps. According to results of our previous work [1], now the rate of tunneling 
across the main energy gap A is either smaller or larger than that for B = 0, depending 
on which magnetic band is initially populated. However, as the first approximation, one 
can neglect this effect and use the standard Landau-Zener equation to estimate the rate of 
tunneling. Numerical results presented in this work undoubtedly confirm this conclusion of 
Ref. (l|. 

The other effect of the magnetic field, which cannot be neglected, is that the tunneled 
particle is accelerated only to a finite kinetic energy, that is in strong contrast with the case 
B = 0. Using classical arguments we estimated the upper boundary for the particle kinetic 
energy, see Eq. f fTT]) . After the acceleration stage, the further dynamics of the classical 
particle is a drift in the direction orthogonal to the vector F, i.e., the Hall current. Our 
numerical simulations of the quantum system clearly indicate the presence of this current. 
Thus, the electric breakdown in the presence of a magnetic field should be interpreted as 
appearance of the Hall current. 
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